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Abstract 

Severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in 
southern China in late 2002 and caused a global outbreak with a fatality rate around 
10% in 2003. Ten years later, a second highly pathogenic human CoV, MERS-CoV, 
emerged in the Middle East and has spread to other countries in Europe, North Africa, 
North America and Asia. As of November 2017, MERS-CoV had infected at least 
2102 people with a fatality rate of about 35% globally, and hence there is an urgent 
need to identify antiviral drugs that are active against MERS-CoV. Here we show that 
a clinically available alcohol-aversive drug, disulfiram, can inhibit the papain-like 
proteases (PL’s) of MERS-CoV and SARS-CoV. Our findings suggest that 
disulfiram acts as an allosteric inhibitor of MERS-CoV PL?” but as a competitive (or 
mixed) inhibitor of SARS-CoV PL’”. The phenomenon of slow-binding inhibition 
and the irrecoverability of enzyme activity after removing unbound disulfiram 
indicate covalent inactivation of SARS-CoV PL?” by disulfiram, while synergistic 
inhibition of MERS-CoV PL?” by disulfiram and 6-thioguanine or mycophenolic acid 
implies the potential for combination treatments using these three clinically available 


drugs. 
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1. Introduction 

Before 2002, human coronaviruses (CoVs) had the reputation of occasionally 
emerging from zoonotic sources and causing mild respiratory tract infections. In late 
2002, however, without any warning, severe acute respiratory syndrome (SARS) 
emerged and spread by coronaviral infection to become a pandemic, mainly in Asia 
but also in other regions, with a fatality rate of 10% (Hilgenfeld and Peiris, 2013). Ten 
years later, when SARS had almost been forgotten, a second highly pathogenic human 
CoV, MERS, caused the severe respiratory syndrome in the Middle East and then 
spreading to other countries due to human activity (Zaki et al., 2012). MERS-CoV has 
infected at least 2100 people with a high mortality rate of 35% since 2012 


(http://www.who.int/csr/don/7-november-2017-mers-saudi-arabia/en/). Because of 


international travel and climate change, we cannot rule out the possibility of the 
emergence of additional highly pathogenic CoVs in the near future (Menachery et al., 
2015; Menachery et al., 2016). Thus, the development of antiviral drugs effective 
against CoVs is urgently needed. 

CoVs are positive-sense single-stranded RNA viruses. After the virion has 
entered the host cell, two polyproteins, ppla and pplab, are directly translated and 
then cleaved by two viral proteases, main protease (M?”) and papain-like protease 
(PL) (Perlman and Netland, 2009). PL’"® is responsible for the cleavage of non- 
structural proteins (nsp) 1, 2 and 3 while MP?” cleaves all junctions downstream of 
nsp4 (Perlman and Netland, 2009). In addition, PL’ can deubiquitinate or deISGylate 
host cell proteins, including interferon factor 3 (IRF3), and inactivate the pathway of 
nuclear factor «-light-chain-enhancer of activated B cells (NF-«B), resulting in the 


immune suppression of host cells (Clementz et al., 2010; Frieman et al., 2009; Yang 
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et al., 2014; Zheng et al., 2008). Due to its multiple roles in viral replication and host 
cell control, PL” is considered a potential antiviral target. 

Disulfiram is a drug which has been approved by the United States Food and 
Drug Administration (FDA) for use in alcohol aversion therapy since 1951 (Bell and 
Smith, 1949; Krampe and Ehrenreich, 2010; Moore et al., 1998). It is known to 
irreversibly inhibit hepatic aldehyde dehydrogenase (Lipsky et al., 2001). Recent 
studies indicate that disulfiram is able to inhibit other enzymes, such as 
methyltransferase, urease and kinase, all by reacting with important cysteine residues, 
suggesting broad-spectrum characteristics (Diaz-Sanchez et al., 2016; Galkin et al., 
2014; Paranjpe et al., 2014). In addition, there has been a clinical trial investigating 
the usage of disulfiram for reactivating latent HIV in order to make it accessible to 
highly active anti-retroviral therapy (Elliott et al., 2015), and the drug has also been 
shown to act as a “zinc ejector” with respect to hepatitis C virus NSSA protein (Lee et 
al., 2016). However, the effect of disulfiram on viral cysteine proteases is still 
unknown. In the present study, we demonstrate that disulfiram is an inhibitor of 
MERS-CoV and SARS-CoV PL”’s, and furthermore that disulfiram acts on MERS- 
CoV and SARS-CoV PL”® via different inhibition mechanisms. Moreover, we 
investigated the synergies between a number of known PL?” inhibitors and 
disulfiram, and our results point to the possibility of using combination treatments 


involving disulfiram and other clinically available drugs against CoVs. 


2. Materials and methods 
2.1. Recombinant protein production — The SARS-CoV PL”® C271A mutation was 


introduced using the QuikChange mutagenesis kit (Stratagene) and was verified by 
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DNA sequencing. The forward primer was 5’- 
gtacactggtaactatcaggcggetcattacactcatata and the reverse primer was 5’- 
tatatgagtgtaatgacccgcctgatagttaccagtgtac. The MERS-CoV and SARS-CoV PL?"’s and 
the SARS-CoV PL?” C271A mutant protein were produced and purified as previously 
described (Chou et al., 2014; Chou et al., 2012; Lin et al., 2014). Briefly, the cultures 
were grown at 37°C for 4 h, then induced with 0.4 mM isopropyl B-D-1- 
thiogalactopyranoside and grown at 20°C for 20 h. The cell pellet was resuspended in 
lysis buffer (20 mM Tris, pH 8.5, 250 mM NaCl, 5% glycerol, 0.2% Triton X-100, 2 
mM f-mercaptoethanol (BME)), lysed by sonication and then centrifuged to remove 
the insoluble pellet. The target protein was purified from the fraction of soluble 
proteins via nickel affinity chromatography, then loaded onto an S-100 gel-filtration 
column (GE Healthcare) equilibrated with running buffer (20 mM Tris, pH 8.5, 100 
mM NaCl, 2 mM dithiothreitol). For the crystallization of SARS-CoV PL?” in 
complex with glycerol, the reductant was removed and 50 uM disulfiram was added 
to each buffer during the purification process. The purity of the fractions collected 
was analyzed by SDS-PAGE and the protein was concentrated to 30 mg/ml using an 


Amicon Ultra-4 30-kDa centrifugal filter (Millipore). 


2.2. Deubiquitination (DUB) assay — The DUB assay was carried out as previously 
described (Cheng et al., 2015; Chou et al., 2008; Lin et al., 2014). The fluorogenic 
substrate Ub-7-amino-4-trifluoro-methylcoumarin (Ub-AFC) (Boston Biochem) was 
added at a concentration of 0.25 uM along with various concentrations of inhibitors 
into 20 mM phosphate (pH 6.5) and each mixture was incubated at 30°C for 3 min. 


After adding 0.2 uM coronaviral PL’, enzymatic activity was determined by 
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continuously monitoring fluorescence intensity at excitation and emission 
wavelengths of 350 and 485 nm, respectively. The data was fitted to obtain ICso 
according to Eq. (1): 

v = vol(1+ICs0'AT)") (1) 

in which v is the initial velocity in the presence of inhibitor at concentration [J] and vo 
is the initial velocity in the absence of inhibitor, while 7 is the Hill constant. 

In addition, to test for the recoverability of activity, coronaviral PL’® was incubated 
with or without 200 uM disulfiram for 1 h and then desalted using a Sephadex G-25 
column. The DUB activity of 0.2 uM treated enzyme was then determined in the 


presence or absence of 5 mM BME. 


2.3. Steady-state kinetic analysis — The peptidyl substrate Dabcyl-FRLKGGAPIKGV- 
Edans was used to measure the proteolytic activity of PL’. Fluorescence intensity 
was monitored at 329 nm (excitation) and 520 nm (emission) and converted to the 
amount of hydrolyzed substrate based on previous studies (Cheng et al., 2015; Chou 
et al., 2008). For inhibition studies, the reaction mixture contained 9-80 uM peptide 
substrate with 0-200 uM disulfiram in 20 mM phosphate (pH 6.5). MERS-CoV PLP? 
at 0.6 uM and wild-type SARS-CoV PL’ and C271A mutant at 0.05 uM was used, 
respectively. After adding the enzyme to the reaction mixture, fluorescence intensity 
was continuously monitored at 30°C. The increase in fluorescence was linear for at 
least 1 min, and thus the slope of the line represented the initial reaction velocity (v). 

The data obtained for the inhibition of MERS-CoV PL” by disulfiram was 
found to best fit a noncompetitive inhibition pattern in accordance with Eq. (2): 


v= KealE][SV/(C + [/Kis) (Ku + [S])) (2) 


while the data obtained for the inhibition of SARS-CoV PL”® by disulfiram was 
found to best fit a competitive inhibition pattern in accordance with Eq. (3) or a mixed 
inhibition pattern in accordance with Eq. (4): 

V = Kea lE][SV/(C1 + [1/Kis) Ku + LS) (3) 

V = kealEJ[SI/((1 + U/Kis) Km + 1 + [/aKis)S]I) (4) 

in which Keat is the rate constant, [E], [S] and [/] denote the enzyme, substrate and 
inhibitor concentrations, and Ky is the Michaelis-Menten constant for the interaction 
between the peptide substrate and the enzyme. Ki; is the slope inhibition constant for 
the enzyme-inhibitor complex and aKis is the slope inhibition constant for the 
enzyme-substrate-inhibitor complex. The program SigmaPlot 12.5 (Systat Software 


Inc., USA) was used for data analysis. 


2.4. Multiple inhibition assay — To characterize the mutual effects of disulfiram and 
other known PL”” inhibitors, the activity of MERS-CoV PL*® was measured with and 
without either 6-thioguanine (6TG) (0 and 15 uM) or mycophenolic acid (MPA) (0 
and 150 uM) in the presence of various concentrations of disulfiram (0-30 uM), and 
that of SARS-CoV PL?’ was measured with and without either 6TG or N- 
ethylmaleimide (NEM) in the presence of various concentrations of disulfiram (0-24 
uM). The concentrations of the peptidyl substrate and MERS-CoV PL”” were 20 and 
0.6 uM, respectively, while those of the substrate and SARS-CoV PL”® were 15 and 
0.05 uM, respectively. Data obtained from the reactions were fitted to Eq. (5): 

v= vo/(1 + [1)/Kj + [J]/K; + W[L/1/0Kik;) (5) 

where v is the initial velocity in the presence of both inhibitors, [/] and [J] are the 


concentrations of the two inhibitors, vo is the velocity in the absence of inhibitors, Kj 
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and Kj are the apparent dissociation constants for the two inhibitors, and a is a 
measurement of the degree of interaction between the two inhibitors (Copeland, 2000; 


Yonetani and Theorell, 1964). 


2.5. Zinc ejection assays — Release of zinc ions from coronaviral PL?’s was 
monitored as the increase in fluorescence emission from the zinc-specific fluorophore 
FluoZin-3 (Thermo Fisher Scientific) (Lee et al., 2016). Briefly, the protein and 
FluoZin-3 were mixed in 20 mM phosphate buffer (pH 6.5) to concentrations of 5 uM 
and 1 uM, respectively, in the presence or absence of 5 uM disulfiram. Fluorescence 
emission was continuously measured at 25 °C using emission and excitation 
wavelengths of 494 nm and 516 nm, respectively, in a PerkinElmer LS50B 


luminescence spectrometer. 


2.6. Thermostability assays — The change in secondary structure of coronaviral PL?'’s 
in the absence and presence of 5 uM disulfiram was continuously measured using 
ellipticity at 222 nm as the temperature was ramped from 30 to 85 °C in a JASCO J- 
810 spectropolarimeter. The protein at 5 uM was dissolved into 20 mM phosphate 


buffer, pH 6.5. The width of the cuvette was 1 mm. 


2.7. Inactivation mechanism — For the inactivation studies, SARS-CoV PL” (0.05 
uM in 20 mM phosphate buffer, pH 6.5) was incubated with different concentrations 
of disulfiram and peptide substrate, and enzymatic activity was traced for 5 min. All 
progress curves recorded showed an exponential course and were analyzed according 


to the following integrated rate equation (Eq. (6)) (Copeland, 2000): 


[P] = vst + [(Vi + Vs)/Kinact] [1 — exp(-Kinactt)] + d (6) 

in which vy; is the initial velocity, vs is the steady-state velocity, and d is the 
displacement on the y-axis. The replot of kinact versus the concentration of disulfiram 
was fitted to a saturation curve according to Eq. (7) (Copeland, 2000): 

Kinact = KmaxL\/(Kinact + [/]) (7) 

in which Kinact is the dissociation constant of the enzyme-disulfiram complex and kmax 


is the maximum inactivation rate constant. 


2.8. Protein crystallization — Crystals of SARS-CoV PL” in complex with BME or 
glycerol were obtained at 22°C by the sitting-drop vapor-diffusion method. For the 
PL”’-BME complex, the protein at 15 mg/ml was incubated with 0.4 mM disulfiram 
for 1 h and then crystallized. Single crystals were grown in reservoir solution 
containing 16% (w/v) PEG 3350 and 0.1 M Bis-Tris propane (pH 8.0). For the PL’’- 
glycerol complex, protein purified with the addition of 50 uM disulfiram into each 
buffer during the purification process was crystallized at 12.5 mg/ml. Single crystals 
were grown in reservoir solution containing 6% (w/v) PEG 8000 and 0.1 M HEPES 
(pH 8.0). All crystals were cryoprotected in reservoir solution supplemented with 
15% and 25% (v/v) glycerol for PL’"°-BME and PL’”-glycerol, respectively, and then 


flash-cooled in liquid nitrogen. 


2.9. Data collection and structure determination — X-ray diffraction data was 
collected at 100 K on the SPXF beamline 15A1 at the National Synchrotron Radiation 
Research Center, Taiwan, ROC using a Rayonix MX300HE CCD detector at a 


wavelength of 1 A. The diffraction images were processed and then scaled with the 
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HKL-2000 package (Otwinowski and Minor, 1997). The structure was solved by the 
molecular-replacement method with Phaser (McCoy et al., 2007) using the structure 
of wild-type SARS-CoV PL’ (PDB entry 2fe8; (Ratia et al., 2006)) as the search 
model. Manual rebuilding of the structure model was performed with Coot (Emsley 
and Cowtan, 2004). Structure refinement was carried out with REFMAC (Murshudov 
et al., 2011). Data-processing and refinement statistics are summarized in Table 3. 
The crystal structures of the SARS-CoV PL?°-BME complex and SARS-CoV PL?®- 
glycerol complex have been deposited in the Protein Data Bank (PDB entries 5y3q 


and 5y3e for PL’°-BME and PL?'°-glycerol, respectively). 


3. Results and discussion 

3.1. The inhibition of MERS-CoV and SARS-CoV PL?’’s by disulfiram — PL?'’s are 
cysteine proteases that use the thiol group of cysteine as a nucleophile to attack the 
carbonyl group of the scissile peptide bond (Chou et al., 2014; Han et al., 2005; 
Verma et al., 2016). Inhibition can be expected if the catalytic cysteine of a PL’ is 
interfered with or modified (Cheng et al., 2015; Chou et al., 2008). Disulfiram is 
known to be a thiol-reactive compound that can covalently modify cysteine residues 
(Diaz-Sanchez et al., 2016; Galkin et al., 2014; Lipsky et al., 2001; Paranjpe et al., 
2014). To determine whether disulfiram can inhibit coronaviral PL’'’s, the DUB 
activity of MERS-CoV and SARS-CoV PL’ was measured in the presence of 
various concentrations of disulfiram. Interestingly, disulfiram showed a dose- 
dependent inhibitory effect on both proteases with ICso values in the micromolar 
range (Fig. 1). Next, to elucidate the kinetic mechanisms of the interactions between 


disulfiram and the two PL?"’s, the proteolytic activity of each enzyme was measured 
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in the presence of various concentrations of a peptidyl substrate and disulfiram. The 
results were then fitted to different kinetic models (competitive, noncompetitive, 
uncompetitive and mixed inhibition). Surprisingly, disulfiram showed a 
noncompetitive inhibition pattern against MERS-CoV PL’® (Fig. 2A) but a 
competitive inhibition pattern against SARS-CoV PL?” (Fig. 2B). This inconsistency 
is quite intriguing since the two enzymes share a similar overall structure and an 
identical catalytic triad (Bailey-Elkin et al., 2014; Chou et al., 2014; Lei et al., 2014; 
Ratia et al., 2006), albeit the inhibition constant (Kj;) of disulfiram for MERS-CoV 
PL” is 4.4-fold higher than that for SARS-CoV PL’® (Table 1). Perhaps this 
discovery should not be surprising given that disulfiram is also a noncompetitive 
inhibitor for Citrullus vulgaris urease with a Kis of 67.6 uM (Diaz-Sanchez et al., 
2016), while its ICs9 for Giardia lamblia carbamate kinase is 0.6-1.4 uM (Chen et al., 
2012). Similarly, a previous study mentions that their compound 4 also has different 
recognition specificity for the two PL’'°s (Lee et al., 2015). Our study once again 
suggests broad-spectrum potency for disulfiram, given the versatility it shows even 


against two coronaviral PLP'’s. 


3.2. Binding synergy analysis of coronaviral PL?” inhibitors — The inconsistent 
inhibitory effect of disulfiram against the two PL”’s suggests that the binding modes 
of disulfiram on the two enzymes may be different. To verify this, multiple inhibition 
assays using disulfiram and other known PL” inhibitors, including 6TG, MPA and 
NEM, were performed (Fig. 3) (Chen et al., 2009; Cheng et al., 2015; Yonetani and 
Theorell, 1964). Interestingly but not surprisingly, we found that disulfiram displays a 


synergistically inhibitory effect with either 6TG or MPA on MERS-CoV PL?", with 
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the lines in the Yonetani-Theorell plots intersecting above the x-axis and a values 
below | in both cases (Fig. 3A and B) (Copeland, 2000). In contrast, in the case of 
SARS-CoV PL”, each of the plots displays two parallel lines and both ao values are 
significantly higher than | (Fig. 3C and D), suggesting that binding of disulfiram and 
of 6TG or NEM are mutually exclusive on SARS-CoV PL? (Copeland, 2000). Since 
6TG is a competitive inhibitor of both PL?'’s (Cheng et al., 2015; Chou et al., 2008), 
the contrasting synergy of disulfiram and 6TG on the two PL”’s confirms the 
inconsistent inhibitory pattern of disulfiram (Figs. 2 and 3). Furthermore, MPA has 
previously been shown to be a noncompetitive inhibitor of MERS-CoV PL?” and to 
work synergistically with 6TG to inhibit MERS-CoV PL? (Cheng et al., 2015). 
Combining those results with our results regarding the binding synergy of disulfiram 
and 6TG or MPA (Fig. 3A and B), we propose that disulfiram may occupy a third 
binding site on MERS-CoV PL”, neither a site at the active center nor the MPA 
binding site. Next, we evaluated PL” inhibition in the presence of disulfiram 
combined with 6TG and/or MPA by proteolytic assays using a peptidyl substrate. We 
found that the ICso of disulfiram against MERS-CoV PL’ showed a 1.6-fold 
decrease in the presence of 15 uM 6TG and a 5.2-fold decrease at 15 uM 6TG when it 
was tested in combination with 150 uM MPA (Table 2). For comparison, in the case 
of disulfiram against SARS-CoV PL", there is no enhanced inhibitory effect in the 
presence of 6TG or NEM. Our results suggest a potential for using the above three 
FDA-approved drugs in combination treatments against MERS-CoV. Incidentally, 
previous studies have suggested that MPA may be used in combination treatments 


with interferon against MERS-CoV (Chan et al., 2013). 
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3.3. Disulfiram may also act as a zinc ejector — Previous studies suggested that 
disulfiram can bind to the zinc-bound cysteines in hepatitis C virus NS5A protein 
(Lee et al., 2016). As there are four cysteines bound to a zinc ion in PL?'’s (Fig. S2C 
and $2D) (Bailey-Elkin et al., 2014; Chou et al., 2014), we performed zinc ejection 
assays to test whether these zinc-bound cysteines may be a candidate for the 
aforementioned “third binding site” occupied by disulfiram on MERS-CoV PL”®. In 
the present study, the zinc-specific fluorophore, FluoZin-3, was used to identify the 
release of zinc ion due to the binding of disulfiram to the enzyme (Fig. 4A). 
Unexpectedly, we observed significant zinc release in the presence of disulfiram not 
only from MERS-CoV PL?” but also from SARS-CoV PL?. This result indicates 
that disulfiram may bind not only to the active site but also to the zinc-binding sites in 
SARS-CoV PL". Following this finding, we tried to fit our inhibitory results to a 
mixed inhibition model (Fig. $1). The two Kj; for the enzyme-substrate and enzyme- 
substrate-inhibitor complexes were 6.0 and 43.8 uM, respectively, showing a 7.3-fold 
difference in the binding affinity for the two putative binding sites (Table 1). This 
significant difference may explain why the inhibitory pattern of disulfiram against 
SARS-CoV PL?” looks more like competitive inhibition. Next, the thermostability of 
the two PL?'’s in the absence and presence of disulfiram was evaluated (Fig. 4B). Not 
surprisingly, the melting temperature of both PL’s decreased 10-15 °C in the 
presence of disulfiram. These results conform to our earlier finding that the release of 


zinc ion can destabilize PL’"° (Chou et al., 2012). 


3.4, Time-dependent inhibition of SARS-CoV PL?”° by disulfiram — Disulfiram is 


known to covalently modify cysteine residues and leave a diethyldithiolcarbamate 
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(DDC) moiety to inactivate the carbamate kinase of Giardia lamblia (Galkin et al., 
2014). In the presence of 5 mM BME, however, the inhibitory effect of disulfiram 
against PL'’s is minor and the ICso is larger than 300 uM (Table 2). This suggests 
that the reductant can protect the enzyme and, therefore, that disulfiram may inhibit 
the enzyme by modifying the cysteine in the catalytic triad (Cys 112-His273-Asp287). 
To further investigate this possibility, the DUB activity of the enzyme was measured 
after incubation with 200 uM disulfiram for 1 h followed by removal of the small 
molecules using a Sephadex G-25 column. This treatment resulted in an 84% loss of 
activity, suggesting irreversible inhibition of SARS-CoV PL?” by disulfiram (Fig. 5A, 
right panel). Similarly, in a previous in vivo study, disulfiram-treated aldehyde 
dehydrogenase showed 77% enzyme inhibition as compared to the activity of the 
control (Lipsky et al., 2001). Next, the disulfiram-treated SARS-CoV PL? was 
incubated with 5 mM BME for 10 min, after which activity was measured to test for 
re-activation. We found that 30% of the enzyme’s activity was restored after 
treatment with BME (Fig. 5A, right panel). The rescuing effect of the reductant 
suggests that the modification was due to the disulfide bonding interaction between 
the enzyme and the inhibitor. However, in the case of MERS-CoV PL’, we found 
that treatment with disulfiram resulted in an irreversible loss of activity which was not 
rescued by the addition of the reductant (Fig. 5A, left panel). Previous studies have 
suggested that the release of zinc ion following treatment with EDTA will lead to a 
62% loss of PL’ activity (Chou et al., 2012). This result is consistent with the effect 
of disulfiram on PL”’s. Also, the inability of the reductant to rescue the DUB activity 
of MERS-CoV PL?”, suggesting that disulfiram cannot influence its active site, is 


compatible with disulfiram’s noncompetitive mode of inhibition of the enzyme. 
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On the other hand, proteolytic assays of SARS-CoV PL”® at various 
concentrations of disulfiram showed dose- and time-dependent decay when enzyme 
activity was measured for 5 min (Fig. 5B). By fitting the data to Eq. 6, different kinact 
values at various concentrations of disulfiram were determined and then plotted 
versus those disulfiram concentrations (Fig. 5C). The saturated curvature suggests a 
slow-binding phenomenon due to covalent inactivation (Copeland, 2000), a 
conclusion supported by the irrecoverability of enzyme activity after disulfiram 
removed (Fig. 5A). Best-fit analysis determined a Kinact of 5.4 UM and a kmax of 0.011 
s! (Fig. 5C and Table 1). Interestingly, the Kina value is close to Kis, indicating that 
disulfiram may inactivate the enzyme very soon after binding. For comparison, 
previous studies have indicated that 6-mercaptopurine and 6TG are also slow-binding 
inhibitors against the same enzyme, albeit enzyme activity was recovered after 


removing the inhibitors (Chou et al., 2008). 


3.5. Proposed binding mechanism of disulfiram to SARS-CoV and MERS-CoV PL?"’s 
— The structure of SARS-CoV PL?” in complex with disulfiram should allow us to 
understand the binding mechanism more clearly. Accordingly, we attempted to 
crystallize SARS-CoV PL”® in the presence of disulfiram. Unfortunately, although 
crystals of the protein were formed in the presence of 0.4 mM disulfiram, the crystal 
structure showed only BME-like electron density near the active-site cysteine with no 
omit electron density shown near the zinc-binding site (Fig. S2A and S2C). BME is a 
reducing agent that is added into the purification buffer to stabilize the protein, and 
which is also known to reverse the effect of disulfiram (Table 2, Fig. 5A and Kitson, 


1975). To avoid this effect, we eliminated all reducing agents from the purification 
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process, added 50 uM disulfiram into all purification buffers, and then attempted to 
crystallize the protein purified under these conditions. Although we were able to grow 
crystals under different crystallization conditions, we again obtained an unexpected 
result, as the only omit electron density near the catalytic site was fitted as a glycerol 
molecule (Figs. S2B and S2D). This result might be due to the crystals having been 
cryoprotected in reservoir solution supplemented with 25% (v/v) glycerol. 
Nevertheless, the binding of BME and glycerol near the active site suggests that the 
active site may be accessible to disulfiram. Next, using the aforementioned two 
complex structures, a disulfiram and a DDC molecule were docked into the glycerol 
and BME binding sites, respectively (Fig. 6). DDC may be able to covalently bind to 
residue Cys112 in a manner similar to that of BME (Fig. 6A), while disulfiram may be 
able to occupy the glycerol site (Fig. 6B). Interestingly, in the docking structure of the 
PL’-disulfiram complex, we can see that one sulfur atom of the disulfide bond of 
disulfiram is within 4 A of residue Cys271 at blocking loop 2 (BL2), which is very 
important for substrate and inhibitor binding (Chou et al., 2014; Ratia et al., 2008). 
For comparison, there is a valine at the same site in MERS-CoV PL?” (Bailey-Elkin 
et al., 2014; Chou et al., 2014). To verify the possible inhibitory effect of disulfiram 
due to binding to residue Cys271, inhibition of the SARS-CoV PL’ C271A mutant 
by disulfiram was measured (Fig. $3). Interestingly, we can see a 4.4-fold increase in 
ICso (Table 2) compared with that for inhibition of wild-type SARS-CoV PL?” by 
disulfiram. In addition, the decrease of the melting temperature of the C271A mutant 
following treatment with disulfiram is 6 °C, lower than that of wild-type SARS-CoV 
PL?” treatment with the same inhibitor (Figs. 4B and 4C). These findings suggest that 


disulfiram may inhibit SARS-CoV PL?” partly via the residue Cys271 and support the 
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reliability of the docking of disulfiram on the glycerol binding site. Based on our 
kinetic and structural results, we propose kinetic mechanism schemes for the 
inhibition of the two PL'’s by disulfiram (Fig. 7). Similar to the mechanism in the 
case of disulfiram-treated urease (Diaz-Sanchez et al., 2016), disulfiram may form a 
covalent adduct with SARS-CoV PL” and then leave a DDC on the active-site 
Cys112, preventing downstream acylation and thereby inactivating the enzyme. In 
contrast, disulfiram shows a noncompetitive inhibitory effect against MERS-CoV 


PL” and can synergistically inhibit that enzyme with 6TG and MPA. 


4. Conclusion 

In this study, we found that disulfiram is, respectively, a noncompetitive and 
competitive (or mixed) inhibitor of MERS-CoV and SARS-CoV PL?’s. Multiple 
inhibition assays also support a kinetic mechanism by which disulfiram together with 
6TG and/or MPA can synergistically inhibit MERS-CoV PL”®, but not, due to its 
competitive mode of inhibition, SARS-CoV PL". On the other hand, the results of 
kinetic assays, continued inactivation after the removal of disulfiram, reactivation by 
reductant, and the phenomenon of slow-binding inhibition suggest that disulfiram 
may act at the active site of SARS-CoV PL”, forming a covalent adduct with residue 
Cys112. Crystal structures of the enzyme in complex with glycerol and BME imply 


that the active site is solvent-exposed and accessible for disulfiram or DDC binding. 
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Figure legends 

Figure 1. Inhibitory effects of disulfiram on coronaviral PL’’s. DUB activity of 
MERS-CoV (A) and SARS-CoV (B) PL” in the presence of disulfiram (6-50 uM) 
was measured. The concentration of fluorogenic substrate (Ub-AFC) was 0.25 uM, 
while the concentration of coronaviral PL’ was 0.2 uM in both cases. The lines show 


best-fit results in accordance with the ICs59 equation (Eq. 1). 


Figure 2. Inhibition of coronaviral PL?’s by disulfiram. The proteolytic activity of 
MERS-CoV (A) and SARS-CoV (B) PL”® were measured in the presence of different 
peptide substrate concentrations (9-80 uM) and various concentrations of disulfiram 
(6-50 uM). The solid lines are best-fit results in accordance with noncompetitive (A) 
or competitive (B) inhibition models. The Rg values are 0.989 and 0.977, 
respectively. The experiments were repeated to ensure reproducibility. Kinetic 


parameters such as Km, Keat and Kis from the best-fit results are shown in Table 1. 


Figure 3. Mutual effects of coronaviral PL’ inhibitors. The activity of MERS- 
CoV PL’ was measured without and with either 6TG (A) or MPA (B) in the 
presence of various concentrations of disulfiram, and that of SARS-CoV PL”® was 
measured without and with either 6TG (C) or NEM (D) in the presence of various 
concentrations of disulfiram. The concentrations of peptidyl substrate and MERS- 
CoV PL’® (A and B) were 20 and 0.6 uM, respectively, while those of peptidyl 
substrate and SARS-CoV PL”® (C and D) were 15 and 0.05 uM, respectively. The 
points are the reciprocals of the initial velocities and the lines are the best fit of the 


data to Eq. 5. The results suggest that the a values for the four experiments (A-D) are 
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0.1, 0.17, 18.2 and 109.3, respectively. 


Figure 4. Effect of zinc ion ejection by disulfiram and its influence on PL?” 
stability. (A) MERS- and SARS-CoV PL”® each was incubated without and with 5 
uM disulfiram. The release of zinc ions from the enzyme was detected as the increase 
of the fluorescence signal of the zinc-specific fluorophore FluoZin-3. (B) and (C) 
Thermostability of MERS-CoV PL”®, SARS-CoV PL’ or SARS-CoV PL?® C271A 
mutant in the absence or presence of 5 uM disulfiram was detected by circular 
dichroism spectrometry. The protein concentration was 0.2 mg/ml. The wavelength 
used was 222 nm and the cuvette pathlength was 1 mm. The right and left dotted lines 
show the melting temperature of SARS-CoV PL” without and with disulfiram, 


respectively. These results indicate that disulfiram destabilized the enzyme. 


Figure 5. Slow-binding inhibition of SARS-CoV PL?” by disulfiram. (A) DUB 
activity of disulfiram-treated MERS- and SARS-CoV PL” in the absence or presence 
of 5 mM B-ME. The enzyme was incubated without or with 200 uM disulfiram for | h 
and the mixture was then desalted using a Sephadex G-25 column. The concentrations 
of fluorogenic substrate (Ub-AFC) and enzyme were 0.25 and 0.2 uM, respectively. 
(B) 0.05 uM SARS-CoV PL? was incubated with different concentrations of 
disulfiram (O uM, closed circles; 2-12 uM, open circles), after which its proteolytic 
activity was measured for 5 min using 15 uM peptidyl substrate. The solid lines are 
best-fit results in accordance with the slow-binding equation (Eq. 6). (C) The 
observed inactivation rate constants (Kinact) from panel B were replotted against 


disulfiram concentration. The solid line is the best-fit result in accordance with the 
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saturation equation (Eq. 7). Kinetic parameters Kinact and Kmax corresponding to the 


best-fit curve are shown in Table 1. 


Figure 6. Binding of disulfiram to SARS-CoV PL”. Overlay of model structure of 
SARS-CoV PL? in complex with DDC (magenta) (A) or disulfiram (orange) (B) 
with the crystal structure of SARS-CoV PL’ in complex with ubiquitin (gray, PDB 
code: 4MOW). DDC and disulfiram are modeled based on the binding sites of BME 
and glycerol, respectively. The red dashed lines show putative polar interactions while 
the black dashed line shows the distance between residue Cys271 and disulfiram as 


4.0 A. 


Figure 7. Schemes of proposed kinetic mechanisms for the inhibition of SARS- 
CoV and MERS-CoV PL?® by disulfiram. The upper diagram denotes enzyme 
catalysis, mixed inhibition and inactivation of SARS-CoV PL?” by disulfiram. The 
lower diagram shows noncompetitive inhibition of MERS-CoV PL” by disulfiram 
and triple inhibition with two other FDA-approved drugs, 6TG and MPA. SH 


symbolizes the thiolate of catalytic triad residue Cys. 
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Table 1. Kinetic parameters of disulfiram inhibition of two coronaviral PL’'’s 


SARS-CoV PL’? 

No inhibitor 195+49* 0.18+0.03° 

Disulfiram 54403 1.1+0.03 

Competitive 18.3423" 0.174001" 46+0.4° 

Mixed inhibition 19.5425" 0.18+0.01° 60+1.1° 
43.8 +5.6° 

C271 mutant 

No inhibitor 2464+3.1° 0.12+0.01° 

MERS-CoV PL”? 

No inhibitor 28.8+4.6° 0.01 +0.0004* 

Disulfiram 30.5+1.8" 0.01+0.0003" 20.1+0.7° 


“ The steady-state kinetic parameters of the PL’'’s were determined according to the 
Michaelis-Menten equation. 

In the presence of disulfiram, the best-fitted kinetic parameters and Kj, were 
determined in accordance with competitive (Eq. 3) or mixed inhibition (Eq. 4) and 
noncompetitive (Eq. 2) inhibition models for SARS-CoV and MERS-CoV PL?”, 
respectively. 

The value is aKj;,, the inhibition constant for the enzyme-substrate-inhibitor 
complex. 


¢ Kinact and Kmax Values are from the best fit to the saturation equation (Eq. 7). 
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Table 2. ICso comparison of disulfiram inhibition of PL 


pro 


of other inhibitors by proteolytic activity assay 


s in the absence or presence 


Enzyme ICs (uM) ICs59 fold decrement 
SARS-CoV PL?” inhibited by 

disulfiram 14.2+0.5 - 
with 6TG (15 uM) 21.8+ 1.0 0.7 
with NEM (4 uM) 18.1 +0.7 0.8 
with BME (5 mM) >300 

SARS-CoV PL?” C271A inhibited by 62.7 + 2.0 - 
disulfiram 

MERS-CoV PL”® inhibited by 

disulfiram 207204" - 
with 6TG (15 uM) 14.5+0.4° 1.6 
with MPA (150 uM) 21.7404 1.0 
with 6TG (10 uM) and MPA (100 uM) 13.7£17.0° By, 
with 6TG (15 uM) and MPA (150 uM) 44+0.2° 5.2 
with BME (5 mM) >300 


** p < 0.05 by Student’s T test. 
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Table 3. X-ray diffraction data collection and refinement statistics 


SARS-CoV PL”°-BME 


SARS-CoV PL?®- 


complex glycerol complex 
Data collection 
Space group C2 C2 
Cell dimensions 
a, b,c (A) 151.4, 33.3, 90.7 151.2, 33.4, 90.9 
a, B, y (°) 90, 125, 90 90, 125, 90 
Resolution® (A) 30-1.65 (1.71-1.65) 30-1.65 (1.71-1.65) 
Rages (%) 4.1 (34.7) 4.7 (45.6) 
I/ol 29.0 (3.6) 26.3 (3.6) 
Completeness (%) 99.7 (98.2) 95.5 (94.8) 
Redundancy 3.6 (3.6) 3.5 (3.7) 
Refinement 
Number of reflections 42,759 (6,082) 41,221 (5,917) 
R factor® (%) 14.7 (16.3) 16.2 (17.7) 
Free R factor“ (%) 18.4 (20.1) 19.9 (21.7) 
Number of atoms 2,994 2,899 
Protein 2,676 2,659 
Ligand/ion 16/6 18/6 
Water 298 216 
B-factors (A’) 
Protein 16.5 27.8 
Ligand/ion 27.0/21.3 34.5/31.8 
Water 28.2 34.8 
rmsd 
Bond length (A) 0.007 0.008 
Bond angles (°) L.3 2 
Ramachandran analysis (%) 
Favored 92.3 93.0 
Allowed Td 7.0 


“ The numbers in parentheses are for the highest-resolution shell. 


b = 
Rows = yy y 
hoi 


1,,-U,) 


Do >, I,; , where I;,; 1s the integrated intensity of a given 


reflection and (J;,) is the mean intensity of multiple corresponding symmetry-related 


reflections. 


“R= DIF ay ob 


/ by F° , where Fp and Fy are the observed and calculated structure 
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ACCEPTED MANUSCRIPT 


factors, respectively. 


“Free R is R calculated using a random 5% of data excluded from the refinement. 
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Highlights: 


Disulfiram, a drug for use in alcohol aversion therapy, can inhibit the papain-like 
proteases of MERS-CoV and SARS-CoV. 

Disulfiram is a noncompetitive inhibitor of MERS-CoV papain-like protease. 
Disulfiram, 6-thioguanine and mycophenolic acid can synergistically inhibit 
MERS-CoV papain-like protease. 

Disulfiram is a competitive inhibitor of SARS-CoV papain-like protease. 
Disulfiram is a slow-binding inhibitor that forms a covalent adduct at the active 
site of SARS-CoV papain-like protease. 


